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We revisit the procedure for comparing the tttt spectral function measured in r decays to that ob- 
tained in e'^e~ annihilation. We re-examine the isospin-breaking corrections using new experimental 
and theoretical input, and find improved agreement between the t~ — > ■K'n'^Vr branching fraction 
measurement and its prediction using the isospin-breaking-corrected e^e^ tt^tt^ spectral func- 
tion, though not resolving all discrepancies. We recompute the lowest order hadronic contributions 
to the muon g — 2 using e^e~ and r data with the new corrections, and find a reduced difference 
between the two evaluations. The new tau-based estimate of the muon magnetic anomaly is found 
to be 1.9 standard deviations lower than the direct measurement. 
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I. INTRODUCTION 



Spectral functions determined from the cross sections 
of e^e~ annihilation to hadrons are fundamental quanti- 
ties describing the production of hadrons from the strong 
interaction vacuum. They are especially useful at low en- 
ergy where perturbative QCD fails to describe the data. 
Spectral functions play a crucial role in calculations of 
hadronic vacuum polarisation (VP) contributions to ob- 
servables such as the effective electromagnetic coupling 
at the mass, and the muon anomalous magnetic mo- 
ment. The latter quantity requires good knowledge of 



the low energy spectral function dominated by the tt 
channel. 

During the last decade, measurements of the tt 
spectral function with percent accuracy became avail- 
able [1-4], superseding older and less precise data. The 
former lack of precision data inspired the search for 
an alternative. It was found [5] in form of accurate 
T~ 7r~7r°i/T- spectral functions [6-9], transferred from 
the charged to the neutral state using isospin symme- 
try. With the increasing e^e~ — *■ tt+tt" experimental 
precision, which today is on a level with the r data, sys- 
tematic discrepancies in shape and normalisation of the 
spectral functions were observed between the two sys- 
tems [10, 11]. It was found that, when computing the 
hadronic VP contribution to the muon magnetic anomaly 
using the r instead of the e+e" data for the 27r and Att 
channels, the observed deviation with the experimental 
value [12] would reduce from 3.3 times the combined ex- 
perimental and estimated theoretical error to less than 
1 [13]. ^ 

In this paper, we include recent r —yn tt^Vt data 
from the Belle experiment [14], and revisit all isospin- 
breaking corrections in this channel taking advantage of 
more accurate data and new theoretical investigations. 



II. TAU DATA 

The r-based evaluation in [10, 11] used the r spec- 
tral functions measured by the ALEPH [7], CLEG [8] 
and OPAL [9] experiments for the dominant hadronic 
decay mode t~ — > i:~'k^Vt. We include here a high- 
statistics measurement of the same decay mode per- 
formed by Belle [14]. Rather different experimental con- 
ditions arc met at the Z centre-of-mass energy (ALEPH, 
GPAL) and at the r(4S') resonance (GLEG, Belle). At 
LEP the T+r^ events can be selected with high efficiency 
(> 90%) and small non-r background (< 1%), thus en- 
suring little bias in the efficiency determination. The 
situation is not as favourable at low energy: because the 
dominant hadronic cross section has a smaller particle 
multiplicity, it is more likely to pollute the r sample and 
strong cuts must be applied, resulting in smaller selection 
efficiency with larger relative uncertainty. Gn the other 
hand, the T{AS) machines outperform LEP in statistics 
for r-pair production: Belle's analysis contains 5.4 mil- 
lion h^TT^VT- candidates (72.2 fb""'^ integrated lu- 
minosity), compared to 81 thousand candidates used by 
ALEPH (including the full LEP statistics accumulated 
on the Z pole). Moreover, GLEG and Belle have an ad- 
vantage for the r final state reconstruction since particles 
are more separated in space. The LEP detectors have to 
cope with coUimated r decay products and the granular- 
ity of the detectors, particularly the calorimeters, plays a 
crucial role. Gne can therefore consider ALEPH/GPAL 
and GLEG/Belle data to be approximately uncorrelated 
as far as experimental procedures are concerned.^ These 



^ Experimental correlations are introduced by common systematic 
errors in the Monte Carlo simulation used. All experiments em- 
ploy the same tau decay and radiative corrections libraries, which 
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FIG. 1: Relative comparison between the r — > tt tt^Ut invariant mass-squared measurements from ALEPH, CLEO, OPAL, 
Belle (data points) and the combined result (shaded band) . 



four data sets arc combined to provide the most precise 
T spectral fmiction, using the newly developed software 
package HVPTools [15]. It transforms the original r data 
and associated statistical and systematic covariance ma- 
trices into fine-grained energy bins (1 MeV). In the com- 
bination, when the value of a bin-wise average exceeds 
the number of degrees of freedom (ndof), the error in the 
averaged bin is rescaled by v^xV^Mof to account for in- 
consistencies, which occur because most experiments are 
dominated by systematic uncertainties. Figure 1 shows 
the relative comparison of the combined data with those 
of each experiment.^ 

The branching fraction fi^jr" for n^Ti'^Vr is ob- 

tained from the measured decay channel h^ir'^Vr 
{^hTz°) by subtracting the non-7r contribution from the 
generic charged hadron mode {h~). The average B^t^o 
value from these experiments and the two other LEP ex- 



periments L3 [17] and DELPHI [18] is (25.847±0.101)%. 
Subtracting from this the current world average value 
(0.428 ±0.015)% for t' K'n^Vr [19], gives B^^o = 
(25.42 ±0.10)%, which is the result used in the following. 



III. ISOSPIN-BREAKING CORRECTIONS 

Historically, the conserved vector current (CVC) rela- 
tion between r and e'^e~ data was considered even be- 
fore the discovery of the r lepton [20, 21]. In the limit of 
isospin invariance, the spectral function of the vector cur- 
rent decay r — *■ X~iyr is related to the e+e^ X'^ cross 
section of the corresponding isovector final state X'^, 



i=i( s 47ra^ 



(1) 



are used for the correction of feed-through from non-h^Tr" final 
states, as well as for the determination of the acceptance and ef- 
ficiency after applying the selection requirements. ALEPH and 
OPAL use however data-driven spectral functions for the feed- 
through corrections, so that the resulting correlations should be 
small. They are hence neglected in the average. 
The total systematic uncertainty from Belle is derived from the 
correlated mass spectra containing revised information with re- 
spect to Table V of Ref. [f4], which has been provided to us by 
Belle [f6]. They differ mostly at the low mass region below the 
p mass peak. 



where s is the centre-of-mass energy-squared or equiva- 
lently the invariant mass-squared of the t final state X, a 
is the electromagnetic fine structure constant, and Vi x~ 
is the non-strange, isospin-one vector spectral function 
given by 
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FIG. 2: Left: Isospin-breaking coiTections from Gem, FSR, /3o(s)//3!l(s) and |i<o(s)/J'-(s)p. Right: Isospin-breaking corrections 
in the ratio of J = 1 components of the form factors \Fo{s)/ F-(s)\'^ due to the vr mass sphtting — m^± — m^o, the p mass 
splitting Snip — nip± — m^o , and the difference STp in the p meson widths. 
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FSR(g) P^js) 
Gem(s) plis) 



F-is) 



(3) 



In Eq. (2), [1 / N x)dN x / ds is the normaUsed invariant 
mass spectrum of the hadronic final state, and Bx- de- 
notes the branching fraction of r — > X~{'^)i't (through- 
out this paper, final state photon radiation is implied 
for T branching fractions). We use for the r mass the 
value mr = (1776.84 ± 0.17) MeV [19], and for the CKM 
matrix element = 0.97418 ± 0.00019 [22], which as- 
sumes CKM unitarity. For the electron branching frac- 
tion we use Be = (17.818 ± 0.032)%, obtained [23] sup- 
posing Icpton universality. Short-distance clcctroweak 
radiative effects lead to the correction S'ew = 1.0235 ± 
0.0003 [10, 24-26]. All the s-dependent isospin-breaking 
(IB) corrections are included in i?iB, and discussed in the 
following for the dominant tttt decay channel. 

The first term in Eq. (3) is the ratio FSR(s)/Gem(s), 
where FSR(s) refers to the final state radiative correc- 
tions [27] in the tt+tt^ channel, and Gem(s) denotes the 
long-distance radiative corrections of order a to the pho- 
ton inclusive t" ■n^'n'^i'r spectrum. Gem('S) includes 
the virtual and real photonic corrections and was calcu- 
lated originally in [28] in the framework of the Resonance 
Chiral Theory [29] . In that work the small axial contribu- 
tions to real photon emission were fixed using the axial 
anomalous terms [30]. A recalculation of Gem(s) was 
presented in [31], where the real photon corrections were 
incorporated via a meson dominance model. Since these 



corrections diverge in the soft-energy limit, a small mass 
must be given to the photon as regularisation. Consis- 
tency however requires that the real photon corrections 
are calculated by summing over all three polarisation 
states of the massive photon [32]. If we include the lon- 
gitudinal polarisation according to Ref. [32], the model- 
independent piece of the radiative corrections changes by 
at most 0.3% close to threshold and rapidly vanishes with 
increasing s. 

The Gem(s) correction used in this analysis is based on 
Ref. [31]. We do not apply, however, any correction for 
the contribution from the square of the 7r(a; — s- 7r°7) am- 
plitude, since it is considered as a background by all ex- 
periments and hence subtracted from the measured spec- 
tral functions. On the other hand, we do keep the inter- 
ference between bremsstrahlung and to amplitudes. The 
resulting Gem(s) function is shown by the solid curve in 
the left-hand plot of Fig. 2. The main numerical differ- 
ence between this correction and that of [28] lies below 
the p peak. Since the origin of the difference is presently 
only partly understood, we assign the full effect as sys- 
tematic uncertainty to the Gem correction. 

The second correction term in Eq. (3), I3q{s) / (3'^{s) , 
arises from the tt^-tt" mass splitting and is important 
only close to the threshold (dotted curve in Fig. 2 (left)). 

The third IB correction term involves the ratio of the 
electromagnetic to weak form factors \Fq{s) / F^{s)\ and 
is the most delicate one. Below 1 GeV, the pion form 
factors are dominated by the p meson resonance, such 
that IB effects mainly stem from the mass and width 



4 



differences between the and p° mesons, and from 
w mixing. The overall effect of this correction is shown 
by the dash-dotted curve in Fig. 2 (left). 

Let us analyse in more detail the IB effects in the form 
factors. A direct calculation of the 27r production am- 
plitudes in e+e^ annihilation and r decays using vector 
meson dominance leads to 



F-{s) 



/p-(5), 



1 



s - im^T^is) 



(4) 
(5) 



where Spi^ is a complex p-Lu mixing parameter. Follow- 
ing [7], two phenomcnological fits to the e'^e" form factor 
data have been performed using the Gounaris-Sakurai 
(GS) [33] and Kiihn-Santamaria (KS) [34] parametri- 
sations"^. For the corresponding mixing strengths and 
phases of the fits we find = (2.00 ± 0.06) x lO^^, 

avgiSf^) = (11.6±1.8)°, and \6f^ \ ^ (1.87±0.06) x lO'^, 
arg((5^) = (13.2±1.7)°, respectively. In both parametri- 
sations, the absorptive parts of the p propagators have 
an explicit energy-dependence of the form —iy^TpO,- (s). 

One of the IB effects is associated with the p meson 
width difference. Within an accuracy of 0.1%, the decay 
widths of the p mesons below y/s ~ 1 GeV arc given 
by their photon inclusive rates into tttt modes [38]. A 
direct calculation of the p 7r7r(7) and ttttj decay rates 
shows that the width difference, STp = FpO — Fp- , is given 
by [38] 



{mps/96TrF^)[l3^{s) - (3^{s)] + (0.45 ± 0.45) MeV, which 
at = 775 MeV gives STp = (-0.61 ± 0.45) MeV. Note 
that if electromagnetic effects were ignored ((5o.- = 0) in 
Eq. (6), we would have ST = —1.06 MeV, which is very 
similar to the cases considered in [5, 28]. 

The second input required to assess the IB effects in 
the form factors is the mass splitting between neutral 
and charged p mesons. Using the expected difference 
nipo - nipo^^^ w 3r{p" e+e-)/(2a) = 1.45 McV, be- 
tween dressed and bare p° mass [39], together with the 
experimental value mp± — nipO = (—0.4 ± 0.9) MeV, ob- 
tained by KLOE from a fit to the cj) tt+tt^tt" Dalitz 
plot [40], one finds S-nip = mp±-mpO = (1.0±0.9) MeV, 
which we use here instead of the degeneracy assumed in 
previous analyses [10, 28]. 

The IB effects in the ratio of / = 1 components of the 
pion form factors (except for p-u mixing) are drawn in 
the right-hand plot of Fig. 2. It is noticeable that the 
effects of photonic corrections and of the 7r^-7r° mass 
difference in the p meson widths largely cancel each other. 

Figure 3 shows the relative difference between the e~^e~ 
and the isospin-breaking-correctcd r spectral functions 
versus s. The relative normalisation is consistent within 
the respective errors and the shape is found in better 
agreement than before [11], despite a remaining devia- 
tion above the p-mass-squared. The discrepancy with 
the KLOE data, although reduced, persists. 
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487r 



[l3l{s){l + So)-l3l{s)il + S^)] , (6) 



where gpT^r is the strong coupling of the isospin-invariant 
pTTTT vertex and (Sq,- denote radiative corrections for 
photon-inclusive p ^ tttt decays, which include p ttttj. 
Contrary to expressions used in previous approaches the 
p meson decay widths in Eq. (6) are independent of the 
photon energy cut-off used to separate the p — > 7r7r(7) 
and p TTTTj rates. In addition to the IB arising 



from the 



mass difference, the radiative correc- 



tions to p — > TTTT and their corresponding radiative rates 
produce a splitting in the p meson widths. For in- 
stance, at y/s = nip = 775 MeV, the width difference 
of Eq. (6) is STp w +0.76 McV, compared to the value 
STp « (-0.42 ± 0.58) MeV used in [5]. The difference 
between the two results is mainly due to the effects from 
radiative corrections (the Jq,- terms in Eq. (6)). Our re- 
sults can also be compared to the one used in [28] , STp = 



^ The fits are performed in the full mass range where the e+e~ data 
are available. This differs from those fits performed in Ref. [35] 
in which the fits were limited to a given single e+e~ experiment 
(with data available only below 1 GeV) with fewer number of 
free parameters. We do not use the Hidden Local Symmetry 
efi^ective model [36] and the effective field theory model [37] as 
these models do not include contributions from the high mass 
resonances such as p' and therefore can only be valid for the 
mass range below about 1 GeV. 



IV. UPDATE OF 



TABLE I: Contributions to a^'^-^'^ [tttt, t] (xlO"^°) from the 
isospin-breaking corrections discussed in Sec. III. Corrections 
shown in two separate columns correspond to the Gounaris- 
Sakurai (GS) and Kiihn-Santamaria (KS) parametrisations, 
respectively. 



Source 



had,LOr 



GS model 



(10-1°) 
KS model 



Sew 
Gem 
FSR 

p-uj interference 

m^± — m^o effect on a 

m^± — m^o effect on Fp 

rripi — TTipO 

7r7r7, electrom. decays 



-12.21 ±0.15 
-1.92 ±0.90 
±4.67 ±0.47 
±2.80 ±0.19 ±2.80 ±0.15 



-7.i 



±4.09 
-5.91 ±0.59 



±4.02 

11+019 
"■11-0.11 

-6.39 ±0.64 



Total 



-16.07 ± 1.22 -16.70 ±1.23 
-16.07 ± 1.85 



The IB corrections applied to the lowest order hadronic 
contribution to the muon g — 2 using r data in the dom- 
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FIG. 3: Relative comparison between e^e~ and r spectral functions, expressed in terms of the difference between neutral 
and charged pion form factors. Isospin-breaking (IB) corrections are applied to r data with its uncertainties, although hardly 
visible, included in the error band. 



TABLE II: The IB-corrected a),^'^'^'^ [■k-k,t] (x10"^°) from the measured mass spectrum by ALEPH, CLEO, OPAL and Belle, 
and the combined spectrum using the corresponding branching fraction values. The results are shown separately in two different 
energy ranges. The first errors are due to the shapes of the mass spectra, which also include a small contribution of 0.11 from 
the r mass and \Vud\- The second errors are due to B.^^o and Be, and the third errors from the isospin-breaking corrections, 
which are partially anti-correlated between the two energy ranges. The last line gives the world average branching fraction and 
also the evaluations of the combined spectra (which are not equivalent to the arithmetic averages of the individual evaluations 
- see text). 



Experiment 


2m^± - 0.36 GeV 


r] (10-") 

0.36- 1.8 GeV 


(%) 


ALEPH 
CLEO 
OPAL 
Belle 


9.46 ± 0.33exp ± 0.058 ± 0.07ib 
9.65 ± 0.42cxp ± 0.17s ± 0.07ib 
11.31 ± 0.76cxp ± 0.15b ± 0.07iB 
9.74 ± 0.28cxp ± 0.15s ± 0.07ib 


499.19 ± 5.20exp ± 2.70b ± 1.87ib 
504.51 ± 5.36exp ± 8.77b ± 1.87ib 
515.56 ± 9.98exp ± 6.95b ± 1.87ib 
503.95 ± 1.90oxp ± 7.84b ± 1.87ib 


25.49 ± O.lOstat ±0.09,yst 
25.44 ± 0.12etat ±0.42,y,t 
25.46 ± 0.17stat ±0.29sy.t 
25.24 ± O.Olstat ± 0.39sy.t 


Combined 


9.76 ± 0.14cxp ± 0.04b ± 0.07ib 


505.46 ± 1.97cxp ± 2.19b ± 1.87ib 


25.42 ± 0.10 



inant tttt channel can be evaluated with 



dNx 



Nxds 



1 



1 



Be 

2s 

2 



ds 



K{s) 



4mJ 



Rmjs) 
Sew 



1 



where K{s) is a QED kernel function [41]. 

The numerical values for the various corrections are 
given in Table I for the energy range between the 2tt 
mass threshold and 1.8 GeV. The present estimate of the 
IB effect from long-distance corrections is smaller than 
the previous one [13, 31], because we now use a G'em('S) 
correction in which the contributions involving the poJTr 
vertex are explicitly excluded (except for its interference 
with the QED amplitude) . Its uncertainty corresponds to 
the difference between the correction used in this analy- 
sis and that from Ref. [28]. The quoted 10% uncertainty 
on the FSR and ttt:^ electromagnetic corrections is an 
estimate of the structure-dependent effects (pion form 
factor) in virtual corrections and of intermediate reso- 



nance contributions to real photon emission [38, 42, 43]. 
The systematic uncertainty assigned to the p-uj interfer- 
ence contribution accounts for the difference in a^,'^'^'^'-' 
between two phenomenological fits, where the mass and 
width of the co resonance are either left free to vary or 
fixed to their world average values. 

Some of the corrections in Table I are parametrisation 
dependent. We choose to take the final corrections from 
the Gounaris-Sakurai parametrisation and assign the full 
difference with respect to the KS results^ as systematic 
error. The total correction for isospin breaking amounts 
to (-16.07 ± 1.85) • 10-1° for a^'"^'^'^[Tnr,T], where all 
systematic errors have been added in quadrature except 
for the GS and KS difference which has been added lin- 
early. This correction is to be compared to the value 



* We do not confirm the significant IB correction difli'erence of the 
KS parametrisation on the p — ut interference with respect to the 
GS parametrisation observed in Ref. [35]. 
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(-13.8 ± 2.4) • 10-1° obtained previously [10]. Since 
the FSR correction was previously included, but not 
counted in the IB corrections, the net change amounts 
to —6.9 X 10-1°, dominated by the electromagnetic de- 
cay correction. 

The corresponding IB-corrected a}lf^'^'^'^[KTT,T\ in the 
dominant tt+tt- channel below 1.8 GcV is given in Ta- 
ble II for ALEPH, CLEO, OPAL, Belle, and for the com- 
bined mass spectrum from these experiments. The eval- 
uation at energy below 0.36 GeV is obtained by fitting an 
expansion in s to the corresponding mass spectrum fol- 
lowing the method introduced in [10]. The comparison 
of the fit with the r data at low energy is shown in Fig. 4 
(left). Good agreement is observed. Indeed, a direct de- 
termination using data gives (10.18 ± 0.98cxp) x 10-^° in 0.958 GeV so that the corresponding aji''^'i^"[7r7r, e+e 



third error stems from the uncertainty in the perturba- 
tive evaluation of the inclusive hadronic cross section in 
the energy ranges 1.8-3.7 GeV and beyond 5 GeV. The 
central value decreases from previously 710.3, obtained 
using incomplete isospin corrections [13] and the super- 
seded combined r spectral function from ALEPH, CLEO 
and OPAL. 

We also re-evaluate the lowest order hadronic contri- 
bution to the muon g ~ 2 using e'^e~ data, updating 
our most recent preliminary result [13] with published 
CMD-2 [2] and KLOE [4] data. The results are given 
in Table III. We have separated the evaluation into 
four distinct energy ranges. The most recent e~^e~ data 
from CMD2, SND and KLOE overlap in the range 0.63- 

had.LOr 



Jcxp 

agreement with the fit-based result of (9.76 ± 0.14oxp) x 
10-1°, which is more precise because of the constraint 
F(0) = 1. The evaluation in the remaining energy re- 
gion is performed directly from a finely (1 McV) binned 
mass spectrum obtained using HVPTools by interpolat- 
ing the original measurements with second order poly- 
nomials (conserving by means of renormalisation the in- 
tegral in each bin before and after interpolation). The 
consistent propagation of all errors is ensured by gener- 
ating large samples of pseudo Monte Carlo experiments. 
The uncertainty due to the interpolation procedure is es- 
timated from a test with a known model to be at most 
0.2 X 10" 1°, which is negligible compared to the other 
systematic uncertainties. It is interesting to compare the 
first and second errors between the experiments. The 
first errors are mainly experimental, but also include 
small contributions from the uncertainties in t mass and 
\Vud\- The second errors are due to Bt^t^o, as measured 
by each experiment, and - to a lesser extent - to Se, for 
which a common value has been used everywhere. Belle 
has the most precise experimental precision on the mea- 
surement of the mass spectrum, whereas ALEPH domi- 
nates the Sttti-o measurement. The result aji^'^''"'-*[7r7r, r] = 
515.2±2.Oexp±2.20±1.9iB (if not stated otherwise, this 
and the following numbers for are given in units of 
10-1°) is obtained from the combined tt~tt^ mass spec- 
trum of ALEPH, CLEO, OPAL and Belle using the world 
average SttttO = (25.42± 0.10)%. This result is consistent 
with the direct average 516.1 ± 1.8oxp ± 2.2b ± 1.9ib; ob- 
tained from the four individual calculations. The ex- 
perimental error from the combined spectrum is slightly 
less precise as it accounts for the incompatibility between 
experiments in certain region of the mass spectrum. 



The contributions to 



^had.LO 



from the tt^^tt 27r° and 



27r+27r- channels below 1.8 GeV arc 21.4 ± 1.3cxp ± 0.6ib 
and 12.3 ± l.Ooxp ± 0.4ib, respectively. This leads to the 
complete r-based lowest order hadronic contribution 

had,LO r J 



705.3 ± 3.9exp ± 0.7,ad ± 0.7qcd ± 2.1ib , 



705.3 ±4.5, 



(7) 



where the second error is due to our treatment of (poten- 
tially) missing radiative corrections in old data included 
in the calculation of the dispersion integral [11], and the 



values can be compared. Agreement is observed between 
CMD2 and SND, while KLOE Hes somewhat lower. To 
account for this, we consider two combinations of the 
e'^e~ data, distinguished by either including or exclud- 
ing the KLOE data. The combination of the data is 
performed also using HVPTools [15], to transform the 
original e^e~ bare cross sections and associated statisti- 
cal and systematic covariance matrices into fine-grained 
energy bins (IMeV), taking into account to our best 
knowledge the correlation within each experiment as well 
as between the experiments. The evaluation in the low 
energy range 2m^±-0.36GeV is performed as for the t 
data by fitting an expansion in s to the combined e~^e~ 
data [10] (right-hand plot of Fig. 4), benefiting from ad- 
ditional space-like precision data [44] . The evaluations in 
the other three energy ranges are obtained by integrating 
directly the combined e+e- cross sections (c/. Table III). 

We find for the difference, between the r 

and e'''e--based evaluations in the dominant tt+tt- chan- 
nel 

11.7 ± 3.5ee ± 3.5r+IB , 
10.6±4.3ee±3.5,+IB, 

where the upper (lower) value is for KLOE data included 
(excluded). The discrepancies amount to 2.4 and 1.9 
times the overall errors, respectively. 

Including the contributions from the other hadronic 
channels [13], we find for the total e"'"e--based lowest 
order hadronic evaluation 



^ had.LOr 



,had,LO 



[e+e-] 



689.8 ± 4.3oxp+rad ± 0.7qcD : 

690.9 ±5. 

■^cxp+rad ± 0.7qcd : 



with total errors of 4.4(5.2) when including (excluding) 
KLOE. Adding the other contributions [13] including the 
latest estimate of the light-by-light scattering (LBLS) 
contribution of 10.5 ± 2.6 [45], we obtain the Standard 
Model predictions (still in 10-i° units) 

a-l^M = 11 659 193.2 ± 4.5 ± 2.6 ± 0.2 , 



[e+e- 



11659177.7±4.4±2.6±0.2, 
11659178. 8±5.2±2.6±0. 2. 
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0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.05 0.1 0.15 0.2 0.25 0.3 0.35 



s (GeV^) s (GeV^) 

FIG. 4: Fit of the pion form factor from 4m^ to 0.3 GeV^ using a third order expansion with the constraints F{0) = 1 and 
using the measured pion charge radius-squared from space-like data [44]. The result of the fit to the r data (left) and to e+e~ 
data (right) is integrated only up to 0.13 GeV"^, beyond which we directly integrate over the data points. 



TABLE III: Evaluated a^'"^'^'~' [7r-7r,e+e ] (xlO ^°) contribution from the e+e data, including and excluding KLOE data. 
The errors correspond to the experimental uncertainties with the statistical and systematical errors added in quadrature (but 



shown separately for individual 


experiments). 






Energy range (GeV) 


Experiment 


Incl. KLOE 


1 (10-"') 

Excl. KLOE 


2m^± - 0.36 


Combined e+e~ (fit) 


9.71±0.12exp 


0.36 - 0.63 


Combined e^e^ 


120.27 ± 1.67cxp 


119.63 ± 1.88exp 


0.63 - 0.958 


CMD2 03 


361.82 ± 2.43stat 


zb 2.10syst 




CMD2 06 


360.17 ± 1.75stat 


i 2.83syst 




SND 06 


360.68 ± 1.38stat 


± 4.67syst 




KLOE 08 


356.82 ± 0.39stat 


zb 3.08syst 




Combined e"'"e~ 


358.51 ± 2.41exp 


360.24 ± 3.02oxp 


0.958 - 1.8 


Combined e'^e^ 


15.02 ± 0.36exp 


15.02 ± 0.39oxp 


Total 


Combined e'^e~ 


503.51 ± 3.47oxp 


504.60 ± 4.33cxp 



where the first errors are due to the lowest order hadronic 
contributions, the second error includes higher hadronic 
orders, dominated by the uncertainty in the LBLS con- 
tribution, and the third error accounts for the uncer- 
tainties in the electromagnetic and weak contributions. 
The predictions deviate from the experimental average. 
a^fP = 11659 208.9(5.4)(3.3) [46, 47], by 15.7 ± 8.2 (r),' 
31.2 ±8.1 (e+e- with KLOE) and 30.1 ±8.6 (e+e" with- 
out KLOE), respectively. 

The lowest order hadronic contribution now reaches an 
uncertainty that is smaller than the measurement error 
and comparable in size with the LBLS uncertainty. Fur- 
ther progress in this field thus requires, apart from con- 
tinuously improved low-energy e^e" cross section mea- 



surements, a more accurate muon .g — 2 measurement and 
LBLS calculation. A compilation of this and other recent 
a^'^ predictions, compared to the experimental value, is 
shown in Fig. 5. 

V. CVC PREDICTION OF B^^o 

The CVC relation (1) allows one to predict the branch- 
ing fraction of a heavy lepton decaying into a G-parity 
even hadronic final state, X~ , via the vector current 
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DEHZ 03 (t) 
-133 ±68 

HMNT 07 (eV) 
-287 ±53 

J 07 (e*e") 

-294 ±65 I- 
This work (x) 

-157 ±52 

This work (eV with KLOE) 
-312±51 >- 

This work (e*e" w/o KLOE) 
-301 ±58 I- 

BNL-E821 (average) 
0±63 



-700 -600 -500 -400 -300 -200 
o oexp 



-100 



100 



x10" 



FIG. 5: Compilation of recently published results for 
(in units of 10~^^), subtracted by the central value of the 
experimental average [46, 47]. The shaded band indicates 
the experimental error. The SM predictions are taken from: 
DEHZ 03 [11], HMNT 07 [48], J 07 [49], and the present r- and 
e+e~-based predictions using r and e'^e~ spectral functions. 



TABLE IV: Contributions to B'^Y'^o ( x 10"^) from the isospin- 
breaking corrections discussed in Sec. III. For those cor- 
rections shown in two separated columns, they correspond 
to the Gounaris-Sakurai and Kiihn-Santamaria parametrisa- 
tions, respectively. 



Source 



AB^vc^ (10-^) 
GS model KS model 



Sew 
Gem 
FSR 

p-uj interference 

m^± — m^o effect on a 

m^± — m^o effect on Fp 

mp± - "ipo^^^ 

7r7r7, electrom. decays 



-h0.57±0.01 
-0.07 ±0.17 
-0.19 ±0.02 
-0.01 ±0.01 -0.02 ±0.01 

±0.19 

-0.22 

±0.08 ±0.08 ±0.09 ±0.08 
±0.34 ±0.03 ±0.37 ±0.04 



Total 



-0.69 ±0.19 ±0.72 ±0.19 
±0.69 ±0.22 



with Sniin being the threshold of the invariant mass- 
squared of the final state in e^e" annihilation. This 
relation was tested ever since the discovery of the r lep- 
ton. In the best known vector channel, the 7r~7r° final 
state, it has attained a precision of better than 1% [13], 
and a discrepancy between B'^Y^o 



of 4.5cr was observed.^ CVC comparisons of r branch- 
ing fractions are of special interest because they are es- 
sentially insensitive to the shape of the r spectral func- 
tion, hence avoiding experimental difficulties, such as the 
mass dependence of the tt" detection efficiency and feed- 
through, and biases from the unfolding of the raw mass 
distribution from acceptance and resolution effects. 



1 1 1 1 1 1 1 1 1 1 1 1 
T decays 






Belie 






' 25.24 ±0.01 +0.39 
CI FO 






' 25.44 ±0.1 2 ±0.42 






_, ALEPH 

£10. 43 ± u. 1 u ± u.uy 






1 DELPHI 






' 25.31 ±0.20 ±0.1 4 






L3 






24.62 ±0.35 ±0.50 






, OPAL 


1 — 




' 25.46 ±0.1 7 ±0.29 






T average 

25.42±0.10 


e*e" CVC 








CMD2 03 










25.03 ±0.29 










CMD2 06 




1 •- 






24.94 ±0.31 










SND 06 










24.90±0.36 










KLOE 08 










24.64±0.29 










e*e" average 

24.78 ±0.28 






1 • 1 









23.5 24 24.5 25 25.5 26 26.5 



27 



B(T"^V^7r"7r") (%) 



FIG. 6: The measured branching fractions for 
■k~-k^Vt [7-9, 14, 17, 18] compared to the predictions from 
the e^e~ -k^-k^ spectral functions, applying the isospin- 
breaking corrections discussed in Sec. HI. For the e^e~ re- 
sults, we have used only the data from the indicated exper- 
iments in 0.63 — 0.958 GeV and the combined e^e~ data in 
the remaining energy domains below rrir . The long and short 
vertical error bands correspond to the r and e^e~ averages 
of (25.42 ± 0.10)% and (24.78 ± 0.28)%, respectively. 

Similar to Aajj'''^''^'-'[7r7r, r], we have evaluated the IB 
corrections to 



2 Tra^m^ j 



X 1 



s 
ml 



2s \ [S'ew 
.Rib 



(9) 



1 



where Smin = ('ti^- ± m^o)^. The results are sum- 
marised in Table IV. The corresponding Sj^^o (Table V) 
is (24.78 ± O.lTcxp ± 0.22ib)% and (24.92 ± 0.21exp ± 
0.22ib)%, based on the combined e~^e~ data, includ- 
ing and excluding the KLOE data, respectively. The 



and Stt-ttO s-t a level 



^ The use of the term standard deviation (cr) in this context re- 
quires caution because the results discussed in this paper are 
mostly dominated by systematic uncertainties with questionable 
statistical properties. 
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TABLE V: Evaluated B'^-'^o (x10~^) from the e^e~ data including and excluding KLOE data, respectively. The errors 
correspond to the experimental uncertainties with the statistical and systematical errors added in quadrature (but shown 
separately for individual experiments). The IB uncertainty of 0.22 is not explicitly quoted for the subcontributions. 



T^lnprcv ran CP ((^pVi 


T^i"vn pri m pttI" 


S^^^o (%) 

Incl. KLOE Excl. KLOE 


m^- + m^o — 0.36 


Combined e+e~ (fit) 


0.03 ± O.OOcxp 


0.36 - 0.63 


Combined e"'"e~ 


1.96 ± 0.03cxp 1.94 ± 0.03exp 


0.63 - 0.958 


CMD2 03 


20.67 ± 0.13etat ±0.12syst 




CMD2 06 


20.58 ± 0.08stat ±0.16By.t 




SND 06 


20.54 ± 0.07stat ± 0.27sy.t 




KLOE 08 


20.26 ± 0.02etat ±0.17syst 




Combined e'^e~ 


20.40 ±0.14cxp 20.56 ±0.17cxp 


0.958 - m-r 


Combined e"'"e~ 


2.39 ± 0.06cxp 


Total 


Combined e^e^ 


24.78 ± 0.17cxp ± 0.22iB 24.92 ± 0.21cxp ± 0.22ib 



first error quoted corresponds to the experimental error 
and the second error due to uncertainties in the isospin- 
breaking corrections. It differs from the r measurement 

by (0.64±0.10^±0.28ee)% and (0.50 ± 0.10^ ± 0.30ee)%, 
respectively, which is still substantial, but less significant 
than the previous result [11. 13]. A graphical compar- 
ison between the IB-corrected S^XP and the measured 
branching fractions ir^ir Vt [7-9, 14] is shown in 

Fig. 6. The B'^'^p results are obtained using the e+e^ 
data from CMD2, SND and KLOE in 0.63-0.958 GeV 
and the combined e~^e~ data in the other energy regions. 

VI. SUMMARY 

We have revisited and updated the isospin-breaking 
corrections to r data in the 2t: mode, incorporating new 
ingredients in the long-distance radiative corrections and 
in the mass and width splittings of mesons that enter the 
pion form factors. We find that the r and e+e~ spec- 



tral functions from CMD-2 and SND are now marginally 
consistent, while a disagreement with the KLOE mea- 
surement remains. The corrected r-based result for the 
Standard Model prediction of the muon g — 2 is now 1.9 
standard deviations lower than the direct measurement, 
coming closer to the e"*"e~ value. Similarly, the predic- 
tion of the T~ — > 7r~7r*'i/T- branching fraction with e~^e~' 
annihilation data exhibits a reduced discrepancy with the 
T measurement. 

We arc indebted to H. Hayashii for providing the correlation 
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bieri, C. Bouchiat, G.V. Fedotovich, E.A. Kuraev, W. Marciano 
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knowledges Conacyt (Mexico) for financial support. This work 
is supported by National Natural Science Foundation of China 
(10491303, 10825524, 10775142), 100 Talents Program of CAS (U- 
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